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C^) \ Abstract. The results are presented of a molecular line survey to search for the spectral signature of infall towards 77 850/jm 

continuum sources believed to be candidate high mass protostellar objects. Up to six different transitions, HCO + J= 1 — > 0, 
J= 3 — » 2 and J= 4 — > 3, ILCO 2i2 - In, N2H + J= 1 — > and H 13 CO + J= 3 — > 2, were observed towards each source. Towards 
the peak of the 850pm emission, N 2 H + was typically strong, with a peak antenna temperature of ~ 1 .5K, with a typical linewidth 
of ~ 2 kms~' . The good agreement between the velocity and velocity width of the N 2 H + and H 13 CO + emission suggests that 
both species are tracing similar material in the sources. With respect to the velocity of the N2H + , there is a statistically significant 
excess of blue asymmetric line profiles in both the HCO + J= 1 — > and H 2 CO transitions. This excess reaches levels similar 
to that seen towards samples of low mass protostars, and suggests that the material around these high mass sources is infalling. 
| We identify 22 promising candidate infall sources which show at least one blue asymmetric line profile and no red asymmetric 

profiles. The infall velocity is estimated to be in the range of 0.1 kms~' to 1 kms~' with an implied mass accretion rate of 
between 2 x 10~ 4 M Q /yr and 10~ 3 M G /yr. 
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<s> 1. Introduction 

c3 ■ 

^ The mass, energy and momentum input from high mass stars (M* > 8M Q ) shapes the physical and chemical properties of the 
• ■ *h . interstellar medium in a galaxy and can control the ISM's evolution. Once formed, the high mass stars, although relatively few 
' in number, rapidly determine a molecular cloud's future evolution, regulating the formation of lower mass stars. Despite their 

^ . importanc e there is as yet no generally a ccepted view of how high mass stars form. Proposed models range from accretion driven 
' formation tMcKee & Tan L l2002l I2003I) . similar to that proposed for the formation of low mass stars, competitive accretion in 
clusters (Bonnell, Vine & Bate 2004) or even the agglomeration of lower mass protostars (Bonnell, Bate & Zinnecker 1998). 

Recent work identifying candidate high mass protostellar objects (e.g Lumsden et al. 2002; Molinari et al. 1996; Sridharan 
et al. 2002) has provided samples of luminous young sources with which to address some of the uncertainties about high mass 
star formation. Identifying and studying the accretion flows which collect the material out of which stars form, either directly or 
indirectly, is an important aspect of understanding high-mass star formation. Evidence for the infall of material around a source 
can be provided by analysis of molecular line profiles. Any ordered motions such as rotation, outflow or infall can produce 
recognisable signatures in the line profiles and their spatial distribution (e.g. Adelson & Leung 1988). However unlike outflow 
or rotation which give rise to both red and blue asymmetric lines, infall produces only blue asymmetric lines. The profiles of 
optically thick lines from infalling material have stronger blue shifted emission than red shifted emission (Leung & Brown 1977; 
Anglada et al. 1987; Zhou 1992; Walker, Narayanan & Boss 1994), provided that the excitation temperature of the molecules 
increases towards the centre of the region. This asymmetry arises as the blue shifted emission from the approaching warm gas on 
the far side of the centre of contraction suffers less extinction than the emission from the receding, nearside, material. Although 
for any one object, outflow or rotation could also produce a blue asymmetric line profile along a particular line of sight to a 
source, for an ensemble of sources the presence of infall should manifest itself as an excess of blue asymmetric line profiles 
compared to red asymmetric profiles. On the other hand, an unbiased sample of sources dominated by outflows or rotation would 
be expected to produce roughly equal numbers of blue and red asymmetric lines. 
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An excess of blue asymmetric line profiles is now well established towards low-mass star forming regions, providing strong 
evidence for infall towards these regions. Not only has it been detected in a variety of different tracers and transitions towards low 
mass protostars but also towards starless cores (Mardones et al. 1997; Lee et al. 1999, 2001, 2004; Williams et al 1999; Park et 
al. 1999; Gregersen & Evans 2000; Gregersen et al. 2000; Di Franceco et al. 2001 ; Belloche et al. 2002). In this paper we present 
the first results of a survey for the infall line asymmetry towards a sample of candidate high mass protostellar objects (HMPOs). 



Table 1. Source Positions. The source numbers are from Williams et al. (2004; WFS) but the tabulated positions are the observed 
positions, which as discussed in the text do not necessary precisely agree with those in WFS. The name of the corresponding 
IRAS point source is also given. 
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The sources observed were originally selected from the IRAS point source catalogue by Srid haran et al. 1 J2002I) based on 
criteria established by Ramesh & Sridharan (1997). The IRAS sources were selected to be bright w ith fluxes of greater tha n 90Jy 
and 500 Jy at 60yum and lOOfim respectively and to have colours in the IRAS bands which match the lWood & Churchwell I Jl989l) 
criteria for UCHII regions but not to have the radio continuum emission expected from a UCHII region. These characteristics 
were choosen to identify sources which are luminous and have dust shells similar to UCHII regions but have not formed an HII 
re gion and so may re present earlier stages in the evolution of mass young stars. These sources have subsequently been studied 
bv lBeuther etail J2002rJall and Williams, Fuller & Sridharan (2004, 2005). 
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Table 2. Observed transitions and telescope parameters. 77 is the telescope beam efficiency, FWHM the full width half maximum 
beam width and Ay the channel spacing o f the correlator. The refe rences for the line frequenci es are 1 - ICaselli et al.ldl995l) . 2 - 
lLovasldl992ll2003h . 3 - online version of lPickett et all Jl998l) . 4 - iGreeersen & Evans I OOOll) 
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Table 3. RMS T* A noise levels in spectra. 
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2. Observations 

To search for evidence of infall, single point observations were made towards the peaks of the 850pm S CUBA sources id entified 
by Williams, Fuller & Sridharan (2004, hereafter WFS) as associated with the IRAS sources selected by Sridha ran et al. 1 (2002). 
The observed positions are listed in Table ^ Since the target positions were derived from an initial reduction of the SCUBA 
images, the positions do not necessarily precisely agree with more accurate positions published in WFS but typically agree 
within a few arcseconds. There are, however, six sources where the positions are discrepant by 14" or more (WFS43, 48, 54, 57, 
77,83). 

The observations were carried out using the 30m telescope of the Institut de Radioastronomie Millimetric (IRAM) 1 near 
Granada, Spain and the 15m James Clerk Maxwell Telescope (JCMT) on Mauna Kea, Hawaii 2 . The frequencies of the lines 
observed, together with other frequency dependent parameters are listed in Table |2] The uncertainties in the line frequencies 
(where known) introduce a maximum uncertainty of ±0.03 km s _I (for the HCO + 1-0) to the derived line velocity. 

1 IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain) 

2 The JCMT is operated by the Joint Astronomy Centre in Hilo, Hawaii on behalf of the parent organizations Particle Physics and Astronomy 
Research Council in the United Kingdom, the National Research Council of Canada and The Netherlands Organization for Scientific Research. 
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At the JCMT the J= 4 -> 3 transition of HCO + and the J= 3 -> 2 transition of H 13 CO + were observed using the RxB3 and 
RxA3 receivers respectively over the periods from April to June 2001 and May to June 2002. The HCO + observations used 
both polarization channels of RxB3 and after checking that the two channels agreed in line intensity and center velocity the 
two channels were averaged together for analysis. A portion of these observations were made using in-band frequency switching. 
This resulted in curved, but smooth baselines which were removed using polynomial fits. The H 13 CO + observations were position 
switched using the single channel of RxA3. The off positions were typically 14' away from the sources at a position (+10', +10') 
from the source. The typical rms noise level reached during the observations was 0.08K for the HCO + and 0. 1 K for the H 13 CO + 
and the typical system temperatures were around 750K and 390K respectively. All the JCMT observations made use of the DAS 
autocorrelator spectrometer. The channel separation used for each of the transitions is listed in Table[2] The JCMT data were all 
initially reduced using the JCMT SPECX package and then imported to the IRAM spectral line reduction package CLASS for 
comparison with the IRAM observations. 



The IRAM observations were carried out between December 2001 and January 2002 and in August 2002. The N2H + and 
HCO + J= 1 — > transitions were simultaneously observed in different polarizations using the two 3mm receivers. The typical 
system temperatures and rms noise levels in the resulting spectra were 120K and 0.07 K, and 130K and 0.07K for the HCO + and 
N2H + respectively. The H2CO line was then observed separately using both the facility 2mm receivers. The two H2CO spectra 
were checked for consistency before being averaged together for analysis resulting in spectra with typical system temperatures 
and noise levels of ~ 190K and ~ 0.07K respectively. During the August 2002 observations the J= 3 — > 2 transitions of HCO + and 
H 13 CO + were observed, one with each of the 1mm receiver systems, simultaneously with H2CO. In no case did the observations 
of H 13 CO + provide a higher signal to noise detection than the JCMT observations and only the JCMT data on this transition are 
presented here and considered in the analysis. The rms noise levels for each of the HCO + and H2CO spectra are given in Table[5] 
All the IRAM observations were made using position switching with the same off positions as used at the JCMT. The sources 
were observed using the autocorrelator spectrometer and during the August 2002 run, the newer VESPA correlator. The IRAM 
data were reduced using the CLASS package. 



Considering both the JCMT and IRAM data, there was only evidence of off-position emission contaminating the observations 
in about four spectra in total and in all cases the emission at the off position was considerably weaker than towards the source, 
easily identifiable and did not affect the determination of the line peak velocity. The pointing and focus were regularly checked 
during the observations with both telescopes. The largest pointing errors seen were ~ 3" and typically the pointing offsets were 
less than 2". 



Spectra towards the sources are shown in Figure 1 and 2. The hyperfine structure of the N2FL transition produces three 
components in the spectra. None of the other transitions have hyperfine structure. 



2.1. Reduction 



The reduction and analysis of the line profiles to look for evidence of infall asymmetry is primarily concerned with the line veloc- 
ities and widths, and we will limit the discussion to just these aspects of the emission. The line intensities and column densities 
which can be infered from the optically thin N2FL and H 13 CO + emission can provide useful information on the circumstellar 
material, but we will delay the discussion of these aspects of the emission until a later work. 
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The N2H + spectra were fitted with the seven hyperfine components given in Table fusing the hyperfine fitting routines in 
CLASS to determine the properties of the emission. The parameters defining the hyperfine components are listed in Tableland 
the results for the fit for each source are given in Table|5] The majority of the spectra are well fitted by a single gaussian velocity 
component, however for some sources, for example WFS2, a single component is not a good model. In Table|5]these sources can 
be identified as the sources for which the residual in the region where the line was fitted is significantly larger than in a region 
of the spectrum free from emission. Most of the spectra are also consistent with optically thin emission, although again there are 
exceptions, for example WFS 13 and WFS 38. 




Fig. 2. Spectra of HCO + J= 3 — > 2 towards the sources. The data is from the IRAM 30m telescope. The dashed vertical lines shows the N2H + 
best fit velocity for each source. For presentation the spectra have been rebinned to 0.22 km/s channels. The intensity axis is in units of T" A in 
Kelvin. 



G. A. Fuller, S. J. Williams and T.K. Sridharan: Infall Towards High Mass Protostellar Objects 



11 



Table 4. Hyperfine structure of N2H + . The velocity offset is measured with respect to the FiF — > F'.F' = 23 —> 12 component 
for which we adopt a frequency of 93.1737767 GHz. The statistical weight of the hyperfine component is given in the column 
labelled Weight. 



Component 


Velocity offset 


Weight 


FiF — > f;f 


(km/s) 




10 -> 11 


6.94 


0.04 


12 12 


5.98 


0.18 


11 10 


5.54 


0.11 


22 -> 11 


0.96 


0.18 


23 12 


0.00 


0.26 


21 11 


-0.61 


0.11 


01 -> 12 


-8.01 


0.11 



The only sources towards which N 2 H + was not detected are WFS 57, WFS 83, and WFS 103. Towards WFS77 the N 2 H + 
line was very weak and only very poorly detected. No H I3 CO + emission was detected towards these sources and these sources 
are not considered any further in the analysis below. The N 2 H + line is also very weak towards WFS 110 and WFS 112 an d these 
sources are also excluded from the anal ysis below. The source WFS 48, whilst observed here and by Willia ms et al. I {2004), 
is not part of the Sridhar an et al. 1 J2002I) ' sample and is not included in any of the analysis discussed below. For consistency in 
probing the material along the line of sight to the central sources in each object, the two remaining objects which were observed 
at (relatively) large offsets from the submillimetre peak (WFS43 and 54; Sec 2) were also excluded from the analysis below. 

The results of fitting single Gaussian components to the H 13 CO + lines for which the line peak intensity is greater than 2.5 
times the rms noise are given in Table|6] Comparing the rms of the residuals of the line fit with the rms noise in an emission free 
region of the spectrum shows that all these spectra are well modelled by these single Gaussian components. The 2.5 times the 
rms noise limit corresponds to a detection in the integrated intensity of more than 10 for a line with a typical width of 2 km s -1 . 

The other transitions show a range of line shapes ranging from Gaussian profiles to profiles with multiple peaks and line 
wings. To measure the line peak velocity of these lines, a parabolic fit over a number of channels around the line centre, covering 
the velocity range out to where the line reached about half of its peak intensity, was used. This method was preferred to using 
Gaussian fits to avoid any line wings influencing the fitted velocities. These parabolic fits were all inspected to confirm that they 
provided a good measure of the velocity of the line peak. The results of these parabolic fits are shown in Table As discussed 
below (§EJ, since the HCO + J= 4 — > 3 transition shows the least range in line shape with the majority of lines being Gaussian, 
these lines were also fitted with single Gaussian components for comparison with the parabolic fits. The results for the line 
velocity from both the parabolic and Gaussian fits are also listed in Table[7] 

3. Velocity Structure 

Since multiply peaked line profiles and line shoulders can be due to the presence of multiple, unrelated velocity components 
along the line of sight, rather than the self absorption expected in the presence of infall, it is important to identify the systemic 
velocity of the gas in which the sources are forming. For this purpose an optically thin spectral line is needed as the line peak 
velocity of such a line will reflect the mass-weighted average velocity along the line of sight. In this study we have two species 
which can serve this purpose: N 2 H + and H 13 CO + . 

For the N 2 H + emission the intensity ratio of the hyperfine components of the transition provide a direct measure of the line 
optical depth. The optical depth derived from the hyperfine fits (Table [5} show that for the majority of sources the line has an 
optical depth of less than 1 . Although there are some sources with larger values, even for these sources it is unlikely that optical 
depth significantly affects the estimate of the systemic velocity as the optical depth at the line peak is at very most 0.4 of this 
total value for broad lines and as little as 0.25 for narrrow lines. Although it is harder to estimate the H 13 CO + optical depth, 
the weakness of the line compared to the HCO + 1-0 emission and its simple Gaussian line profile, consistent in width with the 
N^FFC^ IXTl . suggest the emission is optically thin. 

3.1. Systemic Velocities 

Since we have observations of two tracers w hich diffe rent au thors have used to define source systemic velocities in searches 
for infall dMardones et al. L Il997t iGregersen & Evansl l2000l) . we can compare the properties of these two lines to determine 
whether these species are tracing the same material. One test is to compare the linewidth of the H 13 CO + and that derived from 
the hyperfine fit to the N 2 H + . Figure [3] shows the N 2 H + linewidths and the ratio of the H 13 CO + linewidth to the N 2 H + linewidth. 
The mean linewidth ratio is 1.11 + 0.03 (the quoted uncertainty is the standard error on the mean). Although this result suggests 
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Fig. 3. The left panel shows the distribution of N2H + linewidths of the sources. The vertical dashed line indicates the median 
N2H + linewidth of 2.26 km s _1 . The right panel shows the ratio of H 13 CO + linewidth to N2H + linewidth. The median value of 
1.13 indicated by the dashed vertical line. 
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Fig. 4. Comparison of N2H + and H 13 CO + velocities. The figure shows the difference in velocity between H 13 CO + and N2H + in 
units of the N2H + linewidth. The vertical dashed line marks the median value of -0.02. All the sources for which the H 13 CO + 
was detected with a line peak temperature greater than 2.5 times the rms noise level in the spectrum are included in the figure. 



that on average the linewidth of the H 13 CO + is ~ 10% larger than that of the N2H + , this relatively small difference suggests that 
these two species trace very similar material. Figure|4]shows the difference in line velocity in units of the N2H + line width. The 
mean value of this quantity of -0.03, with a standard error on the mean of 0.02, again showing that in general these lines are 
tracing the same material associated with the forming stars. 

Since these tracers appear to be tracing the same material around the sources, but our observations of N2H + have higher signal 
to noise and we detect it towards a larger fraction of the sources, we choose to adopt the N2H + velocity as the source systemic 
velocity. 



4. Line Asymmetries 

To quantify the blue asymmetry of a line we use the iMardones et al. 1 lll997h asymmetry parameter, 6v, which is defined as the 
difference between the line peak velocity of an optically thick tracer, v(thick), and an optically thin tracer, v(thin), measuring the 
systemic velocity, in units of the optically thin line width, dv(thin), 

v(thick) - v(thin) 

6v = (1) 

dv(thin) 

where we use N2H + as the optically thin tracer. 

Table Q and Figure [5] summarize the asymmetry parameter measurements for the HCO + J= 1 — » 0, J= 3 — > 2, J= 4 — > 3 and 
H2CO spectra, measured with respect to the N2H + velocity and line width. Negative values correspond to blueshifted emission 



Table 5. Hyperfine fits to N2H + spectra. For each source the table gives (in order) the velocity of the emission, the FWHM velocity width of the emission, the line optical 
depth, the line peak temperature in units of T* A , derived from the fit to the hyperfine structure, and the rms in an emission-free region of the spectrum and the rms residual in 
the region of the line after fitting to the hyperfine structure of the line. The values in parenthesis are the one standard deviation uncertainties in the associated quantities. The 
quoted optical depth is the peak optical depth the emission would have in the absence of the hyperfine structure. 



Source 


V 
(km/s) 


dv 

(km/s) 




T 


T peak 

(K) 


Baseline 
rms 
(K) 


Line 
rms 
(K) 


Source 


V 
(km/s) 


dv 
(km/s) 


Tau 


T peak 

(K) 


Baseline 
rms 
(K) 


Line 
rms 
(K) 


2 


-16.10(0.01) 


2.71 


(0.01) 


0.11 


(0.01) 


7.03 


0.07 


0.26 


55 


49.65 ( 0.01) 


3.17 


( 0.02) 


0.10 


0.02) 


2.26 


0.07 


0.08 


5 


0.72 ( 0.03) 


1.00 


( 0.08) 


1.09 


( 1.25) 


0.48 


0.07 


0.07 


59 


77.18 (0.02) 


2.54 


( 0.04) 


0.11 


0.01) 


1.56 


0.07 


0.10 


3 


0.22(0.01) 


1.43 


( 0.05) 


0.47 


( 0.38) 


1.26 


0.07 


0.08 


62 


49.52(0.12) 


1.65 


( 0.28) 


0.10 


0.94) 


0.14 


0.06 


0.05 


6 


6.03 ( 0.09) 


1.61 


( 0.21) 


0.10 


( 0.95) 


0.21 


0.07 


0.07 


64 


10.83 ( 0.02) 


3.14 


( 0.05) 


0.42 


0.19) 


0.93 


0.04 


0.05 


8 


32.56(0.01) 


3.36 


(0.01) 


0.77 


( 0.06) 


3.13 


0.07 


0.21 


65 


85.12(0.01) 


3.04 


( 0.03) 


1.80 


0.14) 


2.29 


0.08 


0.12 


12 


111.03 (0.03) 


2.92 


( 0.08) 


0.35 


(0.31) 


0.95 


0.06 


0.09 


66 


66.04 ( 0.02) 


2.53 


( 0.05) 


1.26 


0.24) 


0.78 


0.04 


0.04 


13 


21.88 (0.01) 


3.53 


(0.01) 


3.41 


( 0.02) 


3.01 


0.06 


0.24 


67 


32.83 ( 0.02) 


4.13 


( 0.04) 


0.30 


0.12) 


1.64 


0.06 


0.07 


14 


33.23 (0.01) 


1.88 


( 0.02) 


0.90 


(0.13) 


0.69 


0.02 


0.02 


68 


55.11 (0.01) 


2.80 


(0.01) 


0.10 


0.01) 


1.08 


0.06 


0.06 


15 


59.96(0.01) 


2.48 


(0.01) 


0.10 


(0.01) 


2.44 


0.06 


0.08 


70 


7.05 ( 0.07) 


2.08 


(0.18) 


0.10 


0.51) 


0.32 


0.07 


0.07 


16 


59.76 (0.01) 


3.09 


(0.01) 


0.11 


(0.01) 


4.00 


0.06 


0.15 


71 


14.05 ( 0.04) 


2.29 


(0.13) 


2.44 


0.78) 


0.45 


0.06 


0.08 


17 


45.12(0.01) 


1.86 


(0.01) 


1.11 


( 0.07) 


3.73 


0.06 


0.13 


72 


3.11 (0.01) 


3.93 


( 0.02) 


0.10 


0.04) 


1.94 


0.06 


0.08 


18 


120.93 ( 0.01) 


2.87 


( 0.02) 


0.10 


( 0.02) 


1.91 


0.06 


0.09 


73 


69.89 ( 0.03) 


3.09 


( 0.04) 


0.73 


0.01) 


1.21 


0.10 


0.10 


19 


43.73 (0.01) 


2.81 


(0.01) 


0.22 


(0.01) 


5.33 


0.06 


0.12 


74 


5.05 (0.01) 


2.93 


( 0.03) 


0.10 


0.05) 


2.04 


0.09 


0.09 


20 


34.85 ( 0.06) 


1.57 


(0.16) 


0.10 


(0.17) 


0.28 


0.06 


0.07 


75 


22.95 ( 0.01) 


1.20 


( 0.02) 


2.83 


0.20) 


2.47 


0.06 


0.10 


21 


34.14 (0.03) 


1.10 


(0.11) 


1.03 


( 1.18) 


0.65 


0.08 


0.08 


79 


22.52 ( 0.01) 


1.80 


(0.01) 


0.91 


0.03) 


6.78 


0.04 


0.18 


22 


84.56(0.01) 


2.45 


( 0.04) 


1.54 


( 0.20) 


2.91 


0.12 


0.15 


80 


29.42 ( 0.01) 


1.74 


( 0.03) 


0.55 


0.22) 


1.50 


0.06 


0.07 


23 


84.25 (0.01) 


2.14 


( 0.03) 


0.68 


(0.16) 


2.84 


0.08 


0.13 


82 


20.42 ( 0.02) 


1.55 


( 0.06) 


4.14 


0.76) 


0.61 


0.04 


0.06 


25 


77.96 (0.01) 


2.59 


(0.01) 


0.10 


(0.01) 


4.71 


0.09 


0.21 


85 


11.21 (0.03) 


1.74 


(0.12) 


0.34 


0.60) 


0.49 


0.05 


0.06 


27 


76.74 (0.01) 


2.78 


( 0.03) 


0.26 


(0.10) 


3.99 


0.09 


0.21 


87 


5.55 (0.01) 


0.82 


( 0.02) 


1.84 


0.52) 


1.11 


0.05 


0.06 


28 


84.29(0.01) 


2.26 


( 0.03) 


1.88 


(0.17) 


2.67 


0.08 


0.18 


88 


5.85 (0.01) 


0.64 


( 0.02) 


2.36 


0.84) 


0.82 


0.05 


0.05 


29 


58.89 ( 0.02) 


2.57 


( 0.04) 


1.96 


( 0.25) 


1.58 


0.08 


0.14 


90 


-3.73 ( 0.01) 


2.02 


( 0.01) 


0.10 


0.01) 


5.06 


0.05 


0.12 


30 


95.51 (0.01) 


3.14 


( 0.03) 


0.80 


( 0.09) 


2.67 


0.07 


0.15 


91 


-2.04 ( 0.01) 


0.95 


( 0.03) 


2.20 


0.74) 


0.62 


0.05 


0.05 


33 


109.70 ( 0.01) 


2.70 


(0.01) 


0.10 


(0.01) 


2.83 


0.07 


0.10 


93 


-1.65 (0.01) 


1.40 


( 0.03) 


0.72 


0.24) 


1.52 


0.06 


0.05 


34 


22.59 ( 0.03) 


2.75 


( 0.09) 


0.59 


( 0.39) 


0.95 


0.08 


0.08 


94 


6.07 (0.01) 


1.98 


(0.01) 


0.36 


0.01) 


5.28 


0.05 


0.14 


35 


25.52 ( 0.09) 


3.46 


( 0.22) 


0.80 


( 0.74) 


0.42 


0.09 


0.07 


95 


8.40 (0.01) 


1.35 


( 0.04) 


0.10 


0.07) 


0.82 


0.05 


0.06 


36 


15.84 ( 0.05) 


1.83 


(0.16) 


2.37 


( 1.10) 


0.48 


0.09 


0.08 


96 


-2.99 ( 0.02) 


1.84 


( 0.06) 


0.59 


0.35) 


1.22 


0.07 


0.08 


37 


105.17 (0.02) 


1.94 


( 0.07) 


1.94 


( 0.50) 


1.14 


0.09 


0.09 


97 


-2.14(0.04) 


2.83 


( 0.08) 


0.10 


0.09) 


0.80 


0.08 


0.08 


39 


96.43 (0.01) 


2.33 


( 0.03) 


0.70 


(0.12) 


1.72 


0.04 


0.09 


100 


11.38 (0.01) 


2.36 


( 0.02) 


1.58 


0.12) 


2.34 


0.06 


0.11 


38 


111.55 (0.02) 


1.40 


( 0.07) 


6.30 


( 1.12) 


0.92 


0.09 


0.13 


99 


10.91 (0.01) 


1.43 


( 0.02) 


1.87 


0.14) 


2.62 


0.06 


0.08 


42 


97.61 (0.01) 


2.32 


( 0.03) 


0.10 


(0.01) 


1.77 


0.07 


0.09 


101 


-18.44 ( 0.01) 


1.56 


( 0.04) 


0.12 


0.01) 


1.44 


0.08 


0.08 


43 


86.86 (0.11) 


2.07 


( 0.34) 


0.10 


( 0.20) 


0.21 


0.07 


0.08 


106 


-47.96 ( 0.01) 


1.92 


( 0.05) 


0.22 


0.26) 


1.45 


0.07 


0.09 


44 


102.60 ( 0.01) 


1.72 


( 0.03) 


0.10 


(0.01) 


2.16 


0.07 


0.12 


107 


-45.70 ( 0.03) 


1.57 


( 0.08) 


1.04 


0.59) 


0.79 


0.07 


0.08 


48 


96.51 (0.02) 


3.00 


( 0.07) 


0.14 


( 0.29) 


0.85 


0.06 


0.07 


108 


-53.08 ( 0.01) 


2.61 


( 0.03) 


0.38 


0.10) 


3.23 


0.08 


0.09 


50 


38.96 ( 0.01) 


2.03 


(0.01) 


0.10 


(0.01) 


2.24 


0.06 


0.08 


109 


-44.41 ( 0.01) 


2.40 


( 0.04) 


0.45 


0.17) 


1.93 


0.08 


0.09 


51 


83.18 (0.01) 


1.68 


( 0.04) 


1.54 


( 0.32) 


1.41 


0.07 


0.09 


110 


-54.45 ( 0.20) 


3.26 


( 0.47) 


0.10 


0.22) 


0.10 


0.05 


0.05 


54 


95.91 (0.04) 


2.56 


(0.10) 


0.10 


( 0.08) 


0.78 


0.09 


0.09 


111 


-17.54(0.03) 


1.22 


( 0.08) 


0.84 


0.82) 


0.61 


0.07 


0.07 


55 


49.65 (0.01) 


3.17 


( 0.02) 


0.10 


( 0.02) 


2.26 


0.07 


0.08 


112 


-18.83 (0.09) 


0.92 


( 0.25) 


0.10 ( 19.42) 


0.16 


0.07 


0.08 


59 


77.18 (0.02) 


2.54 


( 0.04) 


0.11 


(0.01) 


1.56 


0.07 


0.10 





















Table 6. Results from Gaussian fits to H 13 CO + spectra. For each source the table gives (in order) the velocity of the emission, the velocity width of the line, the integrated 
intensity J T* A dv, the line peak temperature in units of T* A from the fit, the rms in an emission-free region of the spectrum and the rms residual in the region of the line after 
fitting the Gaussian. The values in parenthesis are the one standard deviation uncertainties in the associated quantities. 



Source 


V 


dv 


Int. I. 


T 


Baseline 


Line 


Source 


V 


dv 


Int. I 


T 


Baseline 


Line 












rms 


rms 












rms 


rms 




(km/s) 


(km/s) 


(K km/s) 


(K) 


(K) 


(K) 




(km/s) 


(km/s) 


(K km/s) 


(K) 


(K) 


(K) 


2 


-16.09 ( 0.03) 


2.46 ( 0.08) 


2.74 ( 0.08) 


1.05 


0.11 


0.12 


71 


14.06(0.11) 


2.44 ( 0.27) 


0.37 ( 0.03) 


0.14 


0.05 


0.03 


3 


0.13 (0.08) 


1.38 (0.19) 


0.55 ( 0.06) 


0.37 


0.12 


0.12 


72 


2.75 ( 0.09) 


5.16(0.21) 


2.53 ( 0.09) 


0.46 


0.08 


0.07 


6 


5.67 ( 0.05) 


1.60(0.11) 


0.97 ( 0.06) 


0.57 


0.11 


0.10 


75 


22.75 ( 0.04) 


1.75 (0.11) 


1.47 (0.08) 


0.79 


0.12 


0.11 


8 


32.96 ( 0.03) 


4.05 ( 0.07) 


4.91 ( 0.07) 


1.14 


0.08 


0.15 


79 


22.23 ( 0.03) 


2.24 ( 0.09) 


2.79 ( 0.09) 


1.17 


0.13 


0.11 


13 


21.79(0.09) 


3.73 ( 0.20) 


1.58 (0.07) 


0.40 


0.08 


0.08 


80 


29.16(0.09) 


1.96(0.19) 


0.83 ( 0.08) 


0.40 


0.13 


0.11 


14 


33.27 ( 0.02) 


2.13 (0.04) 


3.70 ( 0.05) 


1.63 


0.08 


0.09 


82 


19.99 (0.11) 


1.72 (0.24) 


0.46 ( 0.06) 


0.25 


0.05 


0.04 


16 


59.47 ( 0.05) 


3.73 (0.10) 


3.34 ( 0.08) 


0.84 


0.09 


0.12 


85 


11.08 (0.10) 


2.11 (0.23) 


0.62 ( 0.06) 


0.28 


0.09 


0.07 


17 


45.15 (0.04) 


2.46(0.10) 


1.81 (0.06) 


0.69 


0.08 


0.09 


88 


5.84 ( 0.03) 


0.70 ( 0.06) 


0.29 ( 0.02) 


0.39 


0.06 


0.06 


18 


120.71 (0.10) 


3.03 ( 0.24) 


2.02(0.14) 


0.63 


0.18 


0.19 


87 
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while positive values are redshifted. For comparison, Table0also includes the N2H + and H I3 CO + line velocities as measured by 
hyperfine fitting and Gaussian fitting respectively, together with the velocity of the HCO + J= 4 — » 3 transition also derived from 
a Gaussian fit to the line. Comparing the results from the Gaussian and parabolic fits shows that for the lines which are most 
Gaussian in shape the two methods produce very good agreement in the fitted velocity. Examination of those eight sources where 
the results from the two methods differ by more than 0.3 km s _1 shows that the line profiles are clearly double peaked or skewed. 
The rms difference in velocities over all the sources is 0.26 km s _1 with a mean difference of 0.03 km s _1 , while if only those 65 
sources with a velocity difference of less than 0.5 km s _1 are considered the rms velocity difference drops to 0. 14 km s _1 with a 
mean velocity difference of -0.02 kms" 1 . This good correspondence between the Gaussian and peak fitted velocities, especially 
for sources with nearly Gaussian lines confirms that the peak fitting method produces reliable estimates of the line peak velocity. 

Most mechanisms producing asymmetric line profiles towards sources, e.g. rotation, should produce approximately equal 
numbers of red and blue asymmetric profiles, there being no mechanism to favour one colour of asymmetry over the other. 
On the other hand infall should preferentially produce blue asymmetric profiles (e.g. Anglada et al. 19 87; Zhou 1992; W alker, 
Narayanan & Boss 1994). To quantify whether line profiles of a particular colour dominate in a sample, Mardones et al. 1 dl997l) 
defined a quantity E, the blue excess, 

Metal 

where AWe an d A^red are the number of sources which show blue or red asymmetric lines, respectively, and iVtotai is the total num- 
ber of sources observed. The presence of large blue excesses towa rds samples of low ma ss protostars has provided considerable 
support to the case for the detection of infall around such sources jMardones et al. Lfl997t) . 

To apply the blue excess statistic requires the definition of a range of dv to identify asymmetric sources. Mardonnes et al. 
adopted a criterion of \dv\ > 0.25 to indicate that a line profile was asymmetric and for consistency with their analysis we adopt 
the same criterion. It is important to note however that since the N2H + lines towards this sample of sources is significantly broader 
than towards low mass sources, the choice of this same bound on <5v actually corresponds to a larger absolute velocity shift for 
these high mass objects than for the lower mass sources. Combining the estimated uncertainties on the measured line velocities 
and widths we estimate a typical uncertainty in our measured values of 6\ of 0.02 to 0.05. 

The results of measuring the blue excess for the four transitions which are assumed to be optically thick are listed in Table[S] 
This table shows that for the whole sample of sources observed in HCO + J= 3 — > 2 and J= 4 — > 3 have an equal number of 
blue and red asymmetric lines. However in both the J= 1 — > transiton of HCO + and H2CO, there are considerably more blue 
asymmetric lines than red asymmetric ones, with blue excesses of 15% and 19% for the HCO + and H2CO respectively. In terms 
of the actual number of sources with significant asymmetries for HCO + J= 1 — ♦ there are nearly twice as many blue sources as 
red sources while for H2CO there are over three times as many blue lines as red ones. 

If is of course possible that for any set of observations an excess of blue sources could arise by chance even if the underlying 
distribution of sources has equal number of red and blue asymmetric line profile s. So in order to a ssess the statistical significance 
of the measured blue excesses we have used an exact probability binomial test llConover[IT97Tl) . Given that we have a number 
of sources which have a significant red asymmetry (<5v > 0.25) or blue asymmetry (6v < -0.25), if the sources are randomly 
distributed between red and blue asymmetric, we can calculate the probability that we would observe by chance a blue excess as 
large as actually measured. We do this by calculating the number of ways it is possible to distribute a number of objects equal 
to the total of blue plus red asymmetric profiles into two bins, blue and red, and ask how many of these combinations have a 
number of sources in the blue bin equal to or larger than observed. The ratio of this to the total number of possible combinations 
is then the probabilty that the observed excess arises by chance. This probability is listed in Table [S] as the quantity P. A small 
value of P indicates that it is unlikely that a blue excess as large as observed could arise by chance from an intrinsically uniform 
distribution. 

Table|S]shows that for HCO + J= 1 — > the excess of blue asymmetric sources as compared to red, the binomial test indicates 
only a 6% probability of this asymmetry arising by chance. For the H2CO with a total number of asymmetric sources is 22, 
17 blue and 5 red, the probability is just 0.8% that this excess arises by chance. It appears that for these two transitions, but 
particularly the H2CO, the excess of blue asymmetric lines is statistically significant. 

One point to recall about the sources in this sample is that they span a range of distances from 0.7 to 10.6 kpc ( Srid haran et al. L 

120021) and for some objects there is an ambiguity betw een the source near and far kinematic distances leading to the possibility 
that some sources are at distants as great as 13.6 kpc JSridharanetal.L E002). To limit the possibility of including sources at 
very large distances where the telescope beam may contain emission from material not intimately associated with the high mass 
protostars, we formed a subsample of objects which have confirmed distances of less than 8kpc and have looked at the blue 
excess for these objects alone. The median distance of the sources in this subsample is 3.3 kpc. 

For this distance limited subsample the blue excess seen in HCO + J= 1 — > and H2CO become even more pronounced with 
blue excess of 28% and 23% respectively and probabilities of only 2% and 3% that these arise by chance. Restricting the sample 
even further to include only sources with known distances of less than 6 kpc, giving a mean distance of a source of 2.7kpc, the 
blue excesses become even more significant (Table[8}. In the HCO + J= 1 — » the sample falls to 25 sources, 9 blue and 1 red, 
giving a blue excess of 0.31 and a mean asymmetry of -0.22 + 0.09 while for the PbCO the sample falls to 21 sources, 6 blue 



Table 7. Line velocities and asymmetries. Columns 2, 3 and 4 give the line velocities determined from the hyperfine fit to the N2H + and Gaussian fits to the H 13 CO + and HCO + 
J= 4 — > 3 respectively. The next four columns give the velocities derived from parabolic fits to the lines and the final four columns the line asymmetry parameter derived from 
these velocities. 
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Table 7 continued. 
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Table 8. Blue excess statistics. The table shows for each transition the number of sources with blue asymmetric lines, red 
asymmetric lines and the total number of sources, the blue excess, E, the probability, P, of such an excess and the mean asymmetry 
and the standard error on the mean. 
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Fig. 5. Asymmetry for the four transitions for the sample of objects with known distances of <6 kpc. The dashed lines mark the 
|<5v| < 0.25 region which defines the spectra with no asymmetry. 




and none red, giving an excess of 0.32 with a mean asymmetry of -0.07 + 0.04. The associated probabilities drop to 1% and 2% 
respectively. 

In comparing the statistical significance of the blue excesses in this survey with other surv eys for infall it is imp ortant to 
note that other surveys have not used this exact binomial test which we have used. For example, Mardo nes et al. 1(119971) used the 
Student-t test to compare the distribution of <5v with a normal distrbution with zero mean. This kind of test makes an assumption 
about the nature of the underlaying distribution of asymmetries, including those sources which do not show any asymmetry. 
However there are many reasons why a source may not show an infall asymmetry which do not fit well with the hypothesis 
of smooth distribution of asymmetries assumed by the Student-t test. For example a source might be at a slightly different 
evolutionary stage, or might not quite have the excitation conditions necessary to produce the asymmetry. We therefore prefer 
the binomial test which, once stating a criterion which defines which sources have asymmetric profiles, only considers those 
sources. However, for the sake of comparison applying t he Student t-test to just those sources which have a measured asymmetry 
(which is not how, for example. iMardones et al. I Jl997h used the test) does produce probabilities similar to but somewhat larger 
than those given by our binomial test. For example, for the HCO + J= 1 — > data on the whole sample, the Student t-test gives a 
probability of the observed distribution of blue and red sources arising by chance of 0.08, similar to the 0.06 from the binomial 
test. Applied to the whole sample of asymmetry parameters for the same line in the < 6kpc subsample, the Student t-test gives a 
probability of 0.025, compared to 0.01 for the binomial test applied to only the asymmetric lines. 

As mentioned in Section |2~TI there are a number of sources for which the residuals after fitting the simple single velocity 
component model for the N2H + suggest the presence of additional emission. However the presence of these sources in the sample 
do not significantly distort the measured blue excess. To demonstrate this we have excluded from the sample all the sources for 
which the rms of the single component fit is greater than 1.6 times the rms in the line free region of the spectra and recalculated 
the blue excess statistics. For the remaining sources there are 15 blue and 7 red out of 50 HCO + J= 1 — > spectra giving a blue 
excess of 0.16 with an associated probability of 7%, while for the 47 remaining H2CO spectra, there are 13 blue and 3 red giving 
a blue excess of 0.21 with a probability of 1%. Comparing these values with those in Table[8]shows that the relatively poorly 
fitted sources are having negligible effect on the statistics. The same is true if the 6kpc subsample is examined. 

For both the 8kpc and 6kpc samples, both HCO + J= 4 — » 3 and J= 3 — » 2 show relatively large negative blue excesses, in 
other words, an excess of red asymmetric lines. However in no case does the probability show this excess to be less than a 1 1% 
probability of being due to chance and so these red excesses appear consistent with chance variations due to sampling. 

A comparison of the asymmetry parameter measured towards all the sources observed in both HCO + J= 1 — > and H2CO is 
shown in Figure[6] The asymmetry of these lines is relatively well correlated, but for majority of the sources the asymmetry of the 
H2CO is less than half that of the HCO + perhaps suggesting that these two lines are probing some what different regions around 
the stars. This would be consistent with the different critical densities of these two transitions, 5 x 10 4 cm -3 and 2 x 10 6 cm -3 for 
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the HCO + and H2CO respectively. Interestingly the two transitions which show the least blue asymmetry, HCO + J= 3 — > 2 and 
J= 4 — > 3, also have the highest critical densities (3 x 10 6 cirT 3 and 9 x 10 6 cm -3 respectively) perhaps suggesting stronger infall 
in the lower density, outer regions of these sources. However until the molecular excitation and line formation in the circumstellar 
regions of these sources have been modelled in detail, this should only be regarded as speculation. 



Table 9. Summary of sources with at least one significant line asymmetry. Entries marked R or B indicate red and blue asymmetric 
line profile respectively. A dash (-) indicates a line profile with no significant asymmetry and the absence of an entry indicates 
no observation. Sources marked * are considered the best infall candidates. 



Source 




HCO + 




H 2 CO 


Source 




HCO+ 




H 2 CO 




1 -> 


3^2 


4^3 






1 -> 


3^2 


4^3 




2 






R 




55 * 


B 






B 


5 * 


B 








59 * 


B 


B 




B 


8 


R 


B 


B 


R 


62 


R 








13 


R 




R 


R 


64* 


B 


B 




B 


15 


R 








65 


R 


R 


R 


R 


16 * 


B 




B 


B 


66 * 




B 


B 


B 


17 * 


B 








70 * 






B 
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R 




73 * 


B 


B 


B 


B 


21 * 


B 








79 * 


B 




B 


B 


20 * 


B 






B 


80 






R 




22 * 




B 


B 




82 * 


B 




B 


B 


27 


R 




R 




85 


R 






R 


28 


B 




R 


B 


87 * 


B 








29* 


B 








88 * 


B 




B 


B 


30 


R 








90 






R 




33 


R 




R 


R 


91 






R 




34* 


B 








96 * 


B 






B 


35 


R 








97 




R 






36* 








B 


101 


R 


R 




B 


39* 


B 




B 


B 


109 


B 


R 


R 




51 * 


B 






B 


111 


B 


R 




B 
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Fig. 7. Comparison of the properties of the infall candidates (shaded) an d the other sources in th e sample (outline). The left panel 
shows the distrbution of mass derived from the 850//m observations of Will iams et al. I J2004) for the infall candidates and the 
non-infall sources. The centre panel similiarly compares the distribution o f the IR AS luminosity, assuming the near kinematic 
distances for the sources where the distance is ambiguous JSridharan etal.l Eo02). The right panel shows the distribution of the 
ratio of the N2H + optical depth to the 850/im flux as discussed in Section[5] 



4.1. Comparison with Other Surveys 

The level of blue excess seen in this sample of high-mass objects is comparable to that accepted as evid ence of infall in similar 
surveys of low-mass star forming regions. For example, ove r their whole sample Mardones et al. I dl997l) find an excess of 0.25, 
while for several surveys using HCO + J= 3 -> 2 transitions. Evans l (Eo03) quotes excesses of 0. 30 to 0.31. 

In a recent survey of a different sample of young high mass sources using different tracers IWu & Evans I J2003b found an 
excess of 0.21 to 0.29, again very similar to the values found here. Interestingly there is one object in common between this 
survey and that of Wu & Evans. The source WFS 79 (IRAS 19410+2336) is listed in Wu & Evans as 59.78 + 0.06. While they 
find a velocity for this source of 22.63 + 0.07 km s _1 and an optically linewidth of 2.36 km s" 1 , the N2H + transition is somewhat 
narrower with a width of 1.8 kms^'and a slightly lower velocity of 22.52 + 0.01 kms" 1 . Nevertheless, both surveys find blue 
asymmetric line profiles towards this source with all three transitions observed here significantly blue asymmetric. 



5. Infall Candidates 

The observations of multipe transitions towards the sources in this survey provides a second approach to investigate whether 
there is infall around the sources. Since infall is expected to preferentially produce blue asymmetric lines, it seems reasonable to 
expect that this asymmetry should be present in different transitions observed towards the sources. In particular the presence of 
multiple blue asymmetric transitions and the absence of contradictory information, in the form of red asymmetric profiles, could 
provide a good criterion for identifying sources with infall, independent of the overall excess of blue line profiles in the sample. 

Table|9]summarises those sources which have at least one significantly asymmetric line and from this table we can identify 
those sources we believe to the best infall candidates. We do this by requiring an infall candidate to have at least one blue 
asymm etric transition and no red asymmtric line profiles. This criterion, which is the same as that used by Mar dones et al. I 
( 1997) produces 22 promising infall candidates. The majority of these candidates, 14 sources, in fact have two or more blue 
asymmetric profiles and no red profiles. 

The 22 infall candidates range in luminosity from 10 2 4 L o up to 1O 4 8 L with a median luminosity of 10 L and span a 
range of distances from 0.6kpc to 10.6kpc but only two sources are at distances greater than 5.5 kpc, assuming near kinematic 
distances for sources whose distance ambiguity has not been resolved. At these distances the JCMT and IRAM beams are 
sampling regions 0.02pc to 0.4pc (0.6 to 12xl0 17 cm) in radius around the sources. The submillimetre derived mass of the cores 
around the infall candidates range from several to over a tho usand solar masses wi th a median value of 144M . 

A c omparison of IRAS derived luminosity of the sources JSridharan et al. LI2002I) and the mass derived from the 850/im obser- 
vations (Willi ams etafl 12004 for the infall candidates and the other sources in the sample is shown in FigurefJ] Applying a two 
sample Kolmogorov-Smirnov (K-S) test shows that there is no evidence for a difference in either mass or luminosity between the 
two groups of objects. Similarly K-S tests show there is no significant difference in the N2H + lin ewidths, 850^m colum n density, 
mass to luminosity ratio, 850/im to 450/mi spectral index or cold component dust temperature (Sridharan et al. , 2002) between 
the infall candidates and other sources. Inspection of the radio continuum emission and maser properties ( Sridharan et al. , 2002) 
also shows that the infall candidates do not have any different range of properties to the other sources in the sample. 

However, the K-S test does indicate that there are differences (at a confidence level of about 98%) between the infall candi- 
dates and the remainder of the sample in terms of the optical depth of the N2H + and the 850//m peak flux and, at a somewhat 
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Fig. 8. Line peak intensity ratio R = 7biue/7red as a function of infall velocity i n units of the veloci ty dispersion, V m /<x. The line 
profiles were calculated for five different models using the two layer model of Mvers et al.Nl99d) . 

lower significance level, the 850jum integrated flux. The distributions of these quantities suggest that the infall candidates have 
higher N2H + optical depth and less dust continuum emission than the other sources. Figure0shows the distributions of the ratio 
of the N2H + optical depth to 850yum peak flux, which shows a similiar difference between the infall and non-infall sources with 
the infall candidates having a higher ratio of N2H + to dust emission. Although the statistical significance of this difference is not 
concl usive, if real it has interesting implications about the state of these sources. N2H + is prima rily destroyed by r eacting with 
CO ( iHerbst et al. 1. 1 19 75). leading to high N2H + abundances in regions in which CO is depleted Bergi n et alJ ( 1200 II) . Larger CO 
depletion towards the infall candidates might suggest that these sources are less evolved having not (yet) been as extensively 
heated by their embedded sources as the non-infall sources. A direct test of this will be to measure the CO depletion towards 
these sources. 

5.1. Infall Properties 

If further observations of these candidate sources supports their identification as having infall, then detailed modelling of the 
temperature, density, velocity and indeed chemical structure of the sources will be needed to derive the properties of the in- 
fall. Howeve r an or der of magnitude estimate of the mass accretion rate can be derived from applying the two layer model of 
Mvers et al. 1 dl996l) . In an infalling region the strength of the blue-red asymmetry of a line, as measured by R the ratio of the 
intensity of the blue line peak compared to the red line peak, depends strongly on the infall velocity. Figure[8]shows this intensity 
ratio as a function of the infall velocity (V; n ) in units of the velocity dispersion (<x) calculated using the Myers et al. two layer 
model. 

For lines with small asymmetries, R = T\,\ ue IT K & i 1.5, the infall velocity is relatively tightly constrained. For R = 1.5 the 
infall velocity is between ~ 0.1 and ~ 0.25 times the velocity dispersion. For more asymmetric lines the infall velocity exceeds 
0.25 times the velocity dispersion and can become comparable to, or greater, than the velocity dispersion. For the sources and 
transitions in this survey with blue asymmetric lines which show double peak profiles, the value of R is always in excess of 1 .2 
with values ranging up to 3, implying infall velocities from 0.15 to > 1 times the velocity dispersion. There are also clearly blue 
asymmetric line profiles for which a red peak can not be identifed suggesting even larger infall velocities. Adopting the N2H" 1 " 
linewidth as a measure of the velocity dispersion in the circumstellar material this translates into infall velocities between ~ 0.1 
km s^ 1 and 1 km s _1 with a typical value of ~ 0.2 km s . Adopting a size for the infall region similar to the submillimeter maps 
of WFS, 9 x 10 17 cm and a number density of 5 x 10 4 cirT 3 this range of infall velocities implies mass infall rates of between 
2 x 1 0~ 4 M /yr and 1 0~ 3 M /yr. 

Th is range of infall rates is consistent that values inferred from observations of outflows from these sources jBeuther et all 
2002b), although those rates refer to the region where the stellar winds are generated, much closer the central source that 
the observations presented here. Interestingly the lower end of the range of infall rates is also similar to the rates derived for 
some low-mass sources. F or example, the low lumniosity source NGC 1333 IRAS4A has an infall rate of 1.1 x 10~ 4 M Q /yr 
llDi Francesco et al .112001 1). 

6. Summary 

A survey of up to six transitions of HCO + H 13 CO + H2CO and N2FL has been carried out towards a sample of 77 submillimetre 
sources associated with a well selected sample of IRAS sources believed to be high mass protostellar objects. 

Emission from N2FL was detected towards all the sources except for a few cases where the observation was significantly 
offset from the peak of the submillimetre emission. The strength of the N2FL towards the majority of sources and its weakness 
towards these offset positions suggests that the N2FL and dust continuum emission are tracing similar material. The N2FL is 
mostly optically thin and typically well modelled by a single Gaussian velocity component. The N2FL velocity and linewidth are 
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typically in good agreement with those measured for the same source in H 13 CO + , which is also believed to be optically thin and 
is also well modelled by a single Gaussian component. 

Using the N2FL velocity as a measure of the systemic velocity of the circumstellar material associated with the forming 
stars in these sources, the lines of HCO + and H2CO show a range of profile shapes from near Gaussian, centred at the systemic 
velocity or offset, to double peaked. Analysis of the asymmetry of the line profiles shows that in HCO + J= 1 — > and H2CO 
there is a statistically significant excess of blue asymmetric line profiles compared to red asymmetric profiles. This excess is 
larger if the sample is restricted to exclude the most distant sources. Since the emission from infalling material around a central 
heating source is known to produce systematically blueshifted line profiles, this excess is interpreted as statistical evidence that 
the material around these sources is infalling, as has been argued for similar samples of nearby low mass protostars. 

From the observations of multiple lines, we identify a sample of 22 strong infall candidates none of which show any red 
asymmetric line profiles and have at least one blue asymmetric profile. Confirmation of the presence of infalling material around 
these sources requires further observations, in particular maps of the sources to study the spatial distribution of the blue asymmet- 
ric profiles and higher angular resolution observations to better probe the inner material close to the central sources, together with 
detailed modelling of the sources. Nevertheless, the blue asymmetric line profiles together with density and size estimates from 
submillimetre maps of the sources allow estimates of the range of infall velocity which is found to be between 0.1 kms~' and 1 
kms" 1 and the mass infall rate which is estimated to be between 2 x 10~ 4 M o /yr and 10~ 3 M o /yr . Even the smallest estimate is 
consistent with the mass infall rates needed to form a 1OM in ~ 10 5 years (Mc Kee &Tanll2002l) . 
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